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Abstract 

We investigate the impact of recent limits from LHC searches for supersymmetry and from 
direct and indirect searches for dark matter on global Bayesian inferences of the parameter 
space of the Constrained MSSM. In particular we apply recent exclusion limits from the CMS 
aT analysis of 1.1 fb~^ of integrated luminosity, current direct detection dark matter limit from 
XENONIOO, as well as recent experimental constraints on 7-ray fluxes from dwarf spheroidal 
satellite galaxies of the Milky Way from the FermiLAT satellite, in addition to updated values 
for other non-LHC experimental constraints. We extend the range of scanned parameters to 
include a significant fraction of the focus point/hyperbolic branch region. While we confirm 
earlier conclusions that at present LHC limits provide the strongest constraints on the model's 
parameters, we also flnd that the incidence of the recent exclusion limits from FermiLAT and 
XENONIOO on the posterior pdf is strongly dependent on assumptions about theoretical un- 
certainties. On the other hand, when these uncertainties are not treated in a conservative way, 
the new bounds from indirect detection have the power to significantly constrain the focus 
point/hyperbolic branch region. Their effect is then comparable, if not stronger, to that from 
XENONIOO. We further analyze the effects of one-year projected sensitivities on a neutrino 
flux from the Sun in the 86-string IccCube+DeepCore configuration at the South Pole. We 
show that data on neutrinos from the Sun, expected for the next few months at IceCube and 
DeepCore, has the potential to further constrain the same region of parameter space, and can 
yield additional investigating power for the model. 
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1 Introduction 



The Large Hadron Collider (LHC) at CERN started its run at center of mass energy = 7 TeV 
more than a year ago, and in 2011 published several results based on an integrated luminosity 
of 1.1 fb^"*^ [1, 2j. Several searches were performed, aimed at the discovery of signatures that 
might be an indication of new physics beyond the Standard Model (SM). One of the scenarios 
that are most appealing and well motivated from the theoretical point of view is softly broken 
supersymmetry (SUSY) (see, e.g., Ref. [3]). In the past twenty five years or so there has been 
vast proliferation of models for SUSY-breaking. Among models that are mediated by gravitational 
interactions, the Constrained Minimal Supersymmetric Standard Model (CMSSM) [1] is the most 
popular because of reasonable unification assumptions and a reduced number of parameters. In the 
CMSSM three parameters are defined at the scale of Grand Unification: the universal scalar mass 
mo, the universal gaugino mass ?n.i/2i ^-i^d the universal trilinear coupling Aq. Additionally, tan/3 
(the ratio of the expectation values of the two Higgs doublets), is defined at the electroweak scale, 
while the sign of the Higgs/higgsino parameter /i remains undetermined. 

Recently, there have been many studies, both experimental and phenomenological, of the impact 
of bounds from new physics searches at the LHC on the parameter space (PS) of the CMSSM 
El [71 [8]. In particular, a recent ax analysis \T\ excluded for most MSSM models observation of the 
lightest squarks with mass rrig < 700 GeV, and gluinos with rrig < 580 GeV at the 95% confidence 
level (C.L.). Under sensible assumptions these limits translate into an exclusion bound in the (rn-o, 
mx/2) plane at fixed Aq = 0, tan/3 = 10. The bound is shown in Fig. 8 of Ref. [1]. 

Recently, an updated Bayesian analysis of the CMSSM which included the ax bound from CMS 
was performed by some of the authors of this paper [9j. The analyzed region covered the range 
n^o ^ 2 TeV, mi/2 < 1 TeV, and also included the improved limits on the spin-independent (SI) 
neutralino-proton cross section fip^ from direct detection (DD) of dark matter (DM) with 100.9 
days of five data at XENONIOO [ini HH H]. One of the main findings of [9] was that currently 
LHC limits have the most constraining effect on the PS of the CMSSM, while a further impact of 
the DD limit from XENONIOO is limited after theoretical and astrophysical uncertainties are taken 
into account. 

The study [9] also showed that a combination of the new LHC limits with the previously 
considered usual non-LHC constraints, from particularly flavor physics (5 — t- 57, {g — 2)^, Bs — s- 
fi'^fi^, etc.), dark matter abundance ^x^"^ ^^'^ LEP, results in two broad regions of high posterior 
probability. One is located at large mi/2 and small niQ and agrees reasonably well with the findings 
of studies using a frequentist approach [6]. In addition, a sizable region remains in the focus 
point/hyperbolic branch (FP/HB) region [3l [13], the part of PS at large mo, and mi/2 < rriQ. 

In this paper we look deeper into the FP/HB region and perform a Bayesian analysis of the 
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CMSSM with its soft mass parameters extended to mo < 4 TeV and ^ 2 TeV, which is for each 
parameter twice larger than in [S]. While generally inaccessible to direct LHC searches for SUSY, 
the FP/HB region can be probed by results from DD of DM, and further constraints can be derived 
from recent results on indirect detection (ID). In that respect, in addition to updated constraints 
used in [9], we apply the recent upper bounds on the annihilation cross section of neutralinos 
obtained by the FermiLAT Collaboration |14| from 7-ray fluxes from dwarf spheroidal satellite 
galaxies of the Milky Way (dSphs). Additionally, we discuss the impact of the one- year 95% C.L. 
sensitivities for observation of high energy neutrinos from the Sun at IceCube and DeepCore |15] . 

A full list of constraints applied here and also in [9j will be given below. Recently, the impact on 
the CMSSM parameters of LHC limits and DD (but not ID) limits on DM was also investigated in 
Refs. [Si [71 [8]. Ref. [6] adopted frequentist statistics and [7] performed a fixed-grid scan in contrast 
to our Bayesian one, while Refs. [8] included a different set of constraints. In contrast, the main 
goal of this paper is to analyze the impact of ID searches. 

In Sec. [2] we summarize the main features of Bayesian inference, and describe our scanning 
methodology. We explain our implementation of the LHC and XENONIOO constraints in Sec. |3j In 
Sec. |4]we describe our implementation of the bounds from dSphs into the CMSSM PS, and detail 
our calculation of the IceCube/DeepCore sensitivities. We discuss our results in Sec.jsj and finally 
give our summary and conclusions in Sec. |6j 



2 Methodology 

In Bayesian statistics, for a theory described by some parameters m, experimental observables 
£,{m) can be compared with data d and a posterior probability density function (pdf) p{m\d) can 
be calculated through Bayes' Theorem 

p{d\C{m))TT{m) 

P{m\d) = — , (1) 

where the likelihood p(d|^(m)) gives the probability density for obtaining d from a measurement of 
^, the prior 7r(m) parametrizes assumptions about the theory prior to performing the measurement 
and the evidence p{d) represents the assumptions on the data. As long as one considers only one 
model (as we do in this paper) it is a constant in theory parameters, and thus a normalization 
factor. The Bayesian approach yields a simple and natural procedure for calculating the posterior 
pdf of any limited subset of r variables in PS, ^i=i,..,r C m. One just needs to marginalize, or 
integrate, over the remaining parameters 

p{^i=i,..,r\d) = J p{m\d)d^-'m, (2) 
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where n denotes the dimension of the full PS. To describe our methodology for the Bayesian scan 
we use the same notation as [H]. 

We scan over the CMSSM parameters (mo, ^7ii/2) ^Oj tan /?), which are incorporated in the 
likelihood function according to the procedure detailed in [16]. The Bayesian approach also allows 
a consistent treatment of residual uncertainties in the SM. In the CMSSM, the relevant "nuisance" 
parameters are [16]: the pole mass of the top quark m^°^'^, the mass of the bottom quark at its 
nominal scale in the MS'-scheme, mi,(mi,)^^'^ , and the strong and inverse electromagnetic coupling 
constants at the Z-pole in the M5-scheme, as{Mz)'^'^^ and l/aem{Mz)^'^'^ , respectively. 

We perform our scan in the following ranges: 

100 GeV < mo < 4000 GeV 

100 GeV < mi/2 < 2000 GeV 

-2000 GeV < ^0 < 2000 GeV 

3 < tan /3 < 62 . (3) 

sgn/i is fixed at +1 to accommodate the discrepancy of {g — 2)^ with SM predictions [I7j. The 
nuisance parameters are scanned in the range, 



163.7 GeV < m^"'" < 178.8 GeV 



pole 

3.92 GeV < mi,{mf,)'^^ ^ < 4.48 GeV 



0.1096 < a,(Mz)*^^ < 0.1256 

127.846 <l/aemiMz)^< 127.99. (4) 

We adopt logarithmic priors in mo and mx/2 throughout our analysis, and we refer the reader to the 
relevant literature \16\ [T7] for a discussion on the incidence of linear priors on SUSY PS scans. The 
observables we use to constrain the PS, and their experimental and theoretical errors are reported 
in Table [l| and are the same as given in Ref. |9] . We perform our global scan with the package 
SuperBayeS |18J, which we have modified to our needs. As stated above, the non-LHC constraints 
have been updated with respect to previous papers. To these, we add the LHC constraints, DD 
and ID constraints, which we describe in some detail below. 



3 LHC and DD constraints 

3.1 The ar SUSY search 

In order to calculate the likelihood for the LHC signal we generally follow the procedure detailed in 
Ref. [9]. With respect to the previous paper we extend the efficiency maps to the larger values of 



4 



Observable 


Mean 


Exp. Error Thcor. Error Likelihood Distribution 


Reference 


Non-LHC: 




0.1120 


0.0056 


10% 


Gaussian 


119] 


sin^ 0es 


0.231160 


0.00013 


15.0 X 10~^ 


Gaussian 


m 


Mw 


80.399 


0.023 


0.015 


Gaussian 


m 




30.5 


8.6 


1 


Gaussian 


[iniin] 


Bn{13 X,7) X 10" 


3.6 


0.23 


0.21 


Gaussian 


m 


67^(B„ ^ Tv) X 10'* 


1.66 


0.66 


0.38 


Gaussian 


m 




17.77 


0.12 


2.4 


Gaussian 


m 


Bn{Bs M+M") 


< 1.5 X 10"^ 





14% 


Upper limit — Error fn 




LEP — 95% Limits 


rrih 


> 114.4 





3 


Lower limit — Error fn 


[24 




< /(m^) 








Upper limit — Step fn 


m 




> 50 





5% 


Lower limit — Error fn 


([261I27]) 




> 103.5 (92.4) 





5% 


Lower limit — Error fn 


(|ii[2ll) 




> 100 (73) 





5% 


Lower limit — Error fn 


(HiiiT]) 




> 95 (73) 





5% 


Lower limit — Error fn 


(HSIUZI) 




> 87 (73) 





5% 


Lower limit — Error fn 






> 94 (43) 





5% 


Lower limit — Error fn 


m km) 


LHC CMS ttT 1.1/fb 


analysis 












See text 


See text 





Poisson 




XENONIOO 


of (m^) 


< / (to^) — see text 





1000% 


Upper limit — Error fn 




Nuisance 


l/ae™(Mz)*'^ 


127.916 


0.015 





Gaussian 


m 


pole 


172.9 


1.1 





Gaussian 


m 




4.19 


0.12 





Gaussian 


m 




0.1184 


0.0006 





Gaussian 


[20" 



Table 1: The experimental measurements that constrain the CMSSM's parameters and the Stan- 
dard Model's nuisance parameters. Masses are in GeV. The numbers in parenthesis in the list of 
LEP experimental measurements are weaker experimental bounds, which we use for some sparticle 
mass hierarchies. 
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mo and considered here, and improve the evaluation of the detector efficiency. While in [9] the 
geometrical acceptance only was applied to simulate the CMS detector response, here we apply the 
kinematical cuts detailed in the experimental paper [T] with the fast detector simulator PGS4 [30] , 
whose parameters have been set to the specifications given by the CMS Collaboration. Besides, we 
slightly improve the minimization procedure for the difference in missing transverse energy of the 
pseudo-jets. We have compared our efficiency maps for Aq = and tan/3 = 10 with the results 
of [S] in the overlapping region, finding excellent agreement. 

3.2 XENONIOO bound 

We incorporate in our scan the constraints induced on the SI neutralino-proton cross section £7^^ by 
the 90% C.L. upper bound dp gg which was recently published by the XENONIOO Collaboration [lOl 
[TT\ I12j . As described in the experimental papers, the bound is obtained from the p- value of a 
likelihood function which includes the systematic and statistical error in signal and background, as 
well as the uncertainties on scintillation efficiency and escape velocity; see for details. However, 
the experimental analysis does not consider the astrophysical uncertainties associated with the 
chosen velocity distribution and DM halo profile \31\ [321 [33] , nor the nuclear physics uncertainties 
associated with calculations of the 7r-nucleon sigma-term S,rAr |33], which are in fact dominant. Here 
and below we denote these theoretical uncertainties with r. (In contrast, we label the experimental 
uncertainties with a.) The uncertainties r have been quantified in the literature, and amount up 
to approximately ten times the upper bound on the cross section, r ~ 10 x dp gg |9J . We include r 
in the likelihood, following the procedure described in [16j for the treatment of upper bounds. This 
procedure allows for a systematic treatment of the mismatch between true (^) and approximate (^) 
theoretical quantities; see [16j for details. In the limit of negligible a the likelihood function is given 
by a convolution of a step function, limo— s.op(c, gglO = 0(^p,9o ~ ^^"^ ^ Gaussian function 
parametrizing the theoretical uncertainties, p{C\^, t) =Gauss(.^, r), where in this case ^ = cip^. The 
result is a complementary error function likelihood which smears out the experimental bound at 
the desired level, 

pK%l^,r) = lerfc(^i^^ , (5) 
where the error function erfc(a;) is defined in the standard way. 



erfc(x) = 




Since it is obvious (and has been shown in [9]) that such a large theoretical uncertainty greatly 
reduces the impact of the XENONIOO upper bound, for comparison we consider here the idealized 
case where r is much reduced and is assumed to be comparable in size to c^ gQ. We do this to 
quantify the "bare" effect of the experimental bound on our PS. 
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Two recent papers [8] also considered the impact of the XENONIOO bound on CMSSM PS in the 
Bayesian framework but their treatment of the theoretical uncertainties differs from ours. Unlike 
us, they do not smear out the experimental bound, but instead include the theoretical uncertainties 
in the scan as nuisance parameters and eventually marginalize over them. While in principle this 
is an appropriate procedure in the contest of Bayesian statistics, when uncertainties arise from 
having a number of calculations which are based on different assumptions and methodologies and 
thus yield incompatible results, as is the case with the spread in the SttA^ uncertainty [33], it is 
difficult to believe that they follow any statistical distribution. 

More importantly, as we will see below, even by choosing an extremely optimistic (and, for 
foreseeable future, rather unrealistically) small value of r, say r = fXpgo, the XENONIOO bound 
affects the PS only weakly. This is current LHC limits have much stronger impact on the CMSSM 
parameters [9] . As we will show below, unless one neglected the theoretical uncertainties altogether, 
thus taking the experimental bound at face value as a hard-cut, the impact of the bound cr^ gg on 
the FP region remains limited. 

4 ID constraints 

Besides direct searches at accelerators and underground DM detectors, the third direction in the 
quest for DM discovery is to use detectors to look for uncharacteristic signatures pointing to physics 
beyond the SM in cosmic ray radiation. One method for ID of DM is to exploit searches for 7- 
ray spectra characteristic of DM annihilation \14:\ [35]. [36} [37] . Another one is the search for high 
energy neutrinos from annihilation of weakly interactive massive particles (WIMPs) in the Sun or 
(to minor extent) Earth in neutrino telescopes like IceCube |15j . 

4.1 The FermiLAT 7-ray telescope 

We will focus in this subsection on 7-rays, leaving neutrinos to the next. The FermiLAT detector |35j 
in the Fermi satellite has been used in the past few years on several targets for detection of 7- 
rays from DM annihilation. Both multi- wavelength and line searches have been performed [36] . 
Recently, data was published by the FermiLAT Collaboration which improved significantly the 
previous sensitivities to DM searches. The most luminous source is the Galactic Center (GC) in 
the Milky Way, but it is also subject to higher astrophysical background. Better upper limits on 
the DM annihilation cross section were obtained from the 7-ray fluxes of a few dwarf spheroidal 
satellite galaxies (dSphs) [38l EZl Ej of the Milky Way. They are less luminous but completely 
dominated by DM, with little presence of gas or stars. 

The differential flux of DM annihilation-induced 7-rays from a generic target satellite galaxy is 
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given by 

= f f\2^i^0)]dmdn, (7) 

where N^{E) is the photon energy distribution per annihilation, ((Tf)ann is the velocity-averaged 
pair-annihilation cross section and is the WIMP mass. / is the line-of-sight distance, which 
depends on the opening angle 6 from the center of the galaxy; the integration limits read l± = 
D cos 6 lb \J r"^ — D"^ s\t? 6, where D is the distance to the target galaxy, and rt is the tidal radius 
of the DM halo. p{l) is the DM density distribution in the galaxy, and the angular integration is 
performed about a solid angle Ail calculated by pointing the z-axis of the reference frame in the di- 
rection of the galaxy. The double integral in Eq. ([T]) is commonly referred to as the J-factor, whose 
determination depends on the selected DM galactic profile, and on uncertainties in the stellar veloc- 
ity distribution in the galaxy. In a recent publication [14J, the FermiLAT Collaboration improved 
their previous analysis on dSphs [37J by: (a) including two more dSphs to their original set of eight; 
and (b) incorporating the velocity uncertainties on the J-factor in a likelihood approach. The DM 
density distribution was assumed to follow a cuspy Navarro- Frenk- White (NFW) profile [39j, and 
its normalization and scale factor were obtained by comparisons with stellar velocity observations 
through the likelihood approach. By applying this technique, the Collaboration derived a stronger 
95% C.L. upper bounds on the annihilation cross-sections, (cru) j^gs as a function of m^, for different 
final-state channels (i = bb, fJ-" fJ,~^, t~t^ , W~W^). The bounds are reported in Fig. 2 of Ref. [H]. 

In this paper we investigate the effects of the limits {crv)i^g5 on the CMSSM PS. In order to 
effectively constrain our Bayesian credible regions, we need to be able to invert the procedure that 
the experimental group used to translate an upper limit on the DM-related photon flux into a bound 
on a cross section. The 95% C.L. upper bound on the photon flux as a function of energy ^q^{E) 
was not reported in [14j. The Collaboration obtained it through the energy-binned likelihood 
function given in their Eq. (1). The bound depends non trivially on the flux actually observed by 
the detector, the known background, and the astrophysical uncertainties. 

We can approximately reconstruct <I>95 from any of the given limits {av)i^g^ through Eq. Q, 
since each limit was obtained under the single final state assumption. To illustrate the procedure, 
let us select one final-state channel, say bb. The photon spectrum from bb cascades, N!f^\E), 
does not depend on the nature of the DM particle. Thus, we can eliminate the energy dependence 
by integrating over the FermiLAT spectral sensitivity range, from 0.2 GeV to 100 GeV: Nf^i = 
f N^^\E)dE . Assuming the branching ratio BR^^ = 1, the upper bound on the annihilation cross 
section as a function of is related to the upper bound on the photon flux through Eq. ([7]): 
(cju)^5g5 = 87rm^<I>95/(A'^^5 J), where J denotes the J-factor and $95 = / ^g5{E)dE. 

In the MSSM neutralino pair annihilation proceeds through several final states, both at the tree 
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and loop level. From the above discussion it is easy to see that the 95% C.L. bound on the total 
annihilation cross section (cTu)ann,95 of the MSSM can be obtained from the equality, 

bb,95 

87rm^ 87rm| ' 



from where one obtains 

/ \ \- - / 00 MO " OD /„\ 

(o-l'/ann.gS = ^oQ .... ' (9) 



where we sum the branching ratios BRj = {crv) i / (av) g^^^ over all the 29 channels implemented in 
SuperBayeS. In order to reduce the uncertainties, we performed this procedure three times, using 
the {(7v)ijI (o"w)^-^+^95 and (cu),— ,-+^95 lines. The resulting likelihoods were eventually multiplied 
together. We neglected the VFVF-line because the channel opens up only at rriy- ~ mw, giving less 
constraining power on our PS than the other channels. 

We dealt with the theoretical uncertainties related to halo profiles in two different ways. In |14j 
the J-factor was calculated under the assumption of NFW-like DM density profile. The la- 
uncertainties resulting from the likelihood procedure amount to approximately 3%. The Collabo- 
ration claims that their cuspy profile gives conservative results with respect to smoother, core-like 
profiles. We performed a scan by taking this statement at face value, i.e. we considered the theoret- 
ical uncertainty to be zero, and smeared the 95% C.L. bound by a half Gaussian with uncertainty 
3%. This is the Bayesian equivalent to performing a hard cut on the PS, and results from taking 
the Hmit r in Eq. (3.5) of Ref. [T6] . 

On the other hand, the Collaboration claims that by adopting a core profile for the DM distri- 
bution, the J-factor would either increase or stay the same. The authors of Ref. [40J performed a 
Bayesian scan on a five-dimensional PS in order to constrain different models of DM profile densi- 
ties with observations of velocity dispersion from the known dSphs. They considered a general DM 
profile p (see Eq. Q) parametrized by five parameters, 

7-/3 



p{r) = Ps 



1 + 



r 



(10) 



The five parameters varied in the scan were the density ps calculated at the scale radius r^, and the 
outer, transition and inner slopes, a, (3 and 7 respectively. Their findings show that the posterior 
pdf is very poorly constrained in the space spanned by these five parameters. Thus, the highly cuspy 
NEW profile (q = 1, /3 = 3, 7 = 1), and a flat isothermal core profile (a = 2, /3 = 2, 7 = 0) can fit 
the data equally well. We calculated the J-factor for several different cases, setting the integration 
angle to aint = 0.5°, as provided by [H]. We found that core profiles reduce the J-factor by one or 
two orders of magnitude, in accordance with the findings of [IQ] and |41j . thus rendering the upper 
limits ((Tv)i^95 weaker by comparable amount. We will investigate the importance of this theoretical 
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uncertainty below by including it in the likelihood function employed in our scan in a way similar 
to our treatment of the XENONIOO exclusion bound, as described by Eq. ([s]). One just needs to 
replace ctIqq with {crv)tot,95, and fix the theoretical uncertainty to a value reasonably between one 
and two orders of magnitude. We adopt r = 50 x (o"?;)tot,95- 

4.2 Neutrino Telescopes 

In 2010 the construction of all 86 strings of the IceCube neutrino telescope at the South Pole 
was eventually completed [32] • The instrumented volume covers approximately 1 km^ and deploys 
80 strings uniformly spaced 125 m from one another, which carry 60 optical modules each. The 
energy threshold of the muon effective area is estimated by the experimental Collaboration to be 
approximately 100 GeV [l3]. DeepCore, an additional array of six denser strings is also operational. 
Its performance is parametrized by the muon effective volume Ves- The energy dependence of V^g 
is plotted in [44^, and its parametrization can be found in |45j . The resulting energy threshold is 
~ 35 GeV. 

The IceCube Collaboration published in 2009 the first results on a search for neutrinos from 
neutralino annihilation in the Sun, with 104.3 live days of data-taking in their 22-string configura- 
tion [46j. No excess over the atmospheric background was detected. Recently, the Collaboration 
published a new report with the results from a combined analysis of 812 days of lifetime data from 
AMANDA-II in the years 2001-2006, and 149 days in 2008-09 for the 40-string IceCube detec- 
tor [47]. The data again shows no excess over the background. On the other hand, this implies a 
more stringent upper bound on the spin-dependent (SD) neutralino-proton cross section, and thus 
the data is starting to bite into regions of SUSY PS that were previously probed by DD experiments 
only. While presently available data does not allow to draw any definite conclusions, the results for 
the full configuration (which should be published soon) will have significant constraining power, as 
noted already in Refs. [HI Sll [50]. Moreover, in recent year the study of the detector sensitivity 
has become more elaborate, including increasingly detailed detector simulations [51 ^ 1441 HS] . 

In this paper, we show the impact of the expected one-year 95% C.L. sensitivities for the full 
configuration of IceCube and DeepCore on CMSSM PS, in light of the constraints from the LHC, 
DD, and ID described in the previous sections. With respect to |48| I49j . in addition to a more 
accurate treatment of the most constraining limit from LHC results, we update the effective areas 
and volumes to the most recent values |44l 135] , and investigated the separate impacts of different 
detection channels at IceCube and DeepCore. We briefly explain the formalism below. 

The differential neutrino flux on Earth is given by 
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where Teq is half of the WIMP capture rate, the BRj are the branchmg ratios to the final products 
of neutralino pair annihilation, and Rse is the Sun-Earth distance. Upon arriving on Earth, 
the neutrinos undergo charged-current (CC) interactions in and near the detector volume, thus 
producing a muon flux d^fj_/dE^, which can be detected by the optical modules through Cherenkov 
radiation. In order to realistically approximate the detector response, we follow the procedure 
given in Appendix B of Ref. [M] . The muon events observed by the 80-string IceCube detector can 
either be labelled as upward (those that come from neutrinos undergoing CC-scattering outside 
the 1-km'^ instrumented volume), or as contained (those for which the interaction happens inside 
the detector). We denote the differential flux of upward muons detected by the optical modules 
with d^fi/dEj^, where the index "f is added to distinguish it from the muon flux at the interaction 
point, d^fj,/dE^. The difference arises because of energy losses due to ionization, bremsstrahlung, 
pair production and photo-nuclear effects [52] that affect the spectrum of muons transiting in ice. 
The upward event rates at the detector are given by 

= I ' (12) 

where Aef[{Eji) is the muon effective area parametrized in Ref. ^43j. We decided to neglect the 
angular dependence of the effective area since, for the signal, it integrates to a factor approximately 
equal to one. The expected signal rate has to be compared with the atmospheric background, which 
we adopted from Ref. N^q ~ 6.1 yr~^. 

The contained rate, by definition, does not suffer from propagation related losses, and is given 

by 

^events = / ^ ( 1 km^ ) di?^ , (13) 

where the threshold E'thr is conservatively fixed at 100 GeV. The background is ~ 15.6 yr^^. 

Notice that the upward and contained event rates and atmospheric backgrounds are calculated 
for the time of the year the Sun spends below the horizon, to reduce the background from down 
going muons. So, the rates are calculated for six months time, between the March and September 
Equinoxes [HI US]. For DeepCore, the larger instrumented volume outside the detector will be used 
to veto atmospheric muons from above the horizon, so the event rates are calculated over a twelve 
month period. 

The events rate at DeepCore is given by 

A^e°S.ts = n ^^V,s{E,)dE, , (14) 

where the typical scale of the detector is L = 1 km, and Vcf[{E^) is the effective volume shown 
in [H]. The background is N^q ~ 2.5 yr~^. 
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We calculated the signal using a modified version of SuperBayeS, where we updated the muon 
effective area and effective volume to the values cited above. Notice that SuperBayeS (following 
DarkSusy [53]) uses tables from WimpSim [5l] for the propagation of the neutrino flux through 
the Sun and to the Earth, so that matter effects and neutrino oscillations are properly accounted 
for. Also, the signal is calculated without the approximation of thermodynamical equilibrium [53], 
so that the rates slightly depend on ((Tf)ann- We fixed the local DM density to 0.3 GeV/cm^, 
and the root-mean-square of the velocity dispersion in the halo to 270 km/s. We compared our 
spectra and rates with the code used for the calculations performed in [H], and found very good 
agreement. Our 95% C.L. (or 2a) projected sensitivities are calculated in the standard way, with 

-/Vevents/ V^BG = 2. 

5 Results in the CMSSM 

In this section we will present our numerical results. First, to set the ground, we will display regions 
of the CMSSM mass parameters favored by current data from the usual non-LHC constraints, from 
the ot constraint at the LHC, as well as from the XENONIOO limit for which we will consider two 
limiting cases of the theoretical error. Next, we will add the constraint from dwarf spheroidals and 
again will discuss the impact of assuming different theoretical uncertainties. Finally we will discuss 
the projected sensitivities at IceCube/DeepCore. 

In all plots we denote as "non-LHC" the usual set of relevant constraints as summarized in 
the first part of Table [T| while the ax bound is indicated with "LHC (ot)"- The constraints from 
DD and ID are labeled according to the size of the theoretical uncertainties assumed for them (see 
Sees. [3] and [4]). In particular, for the XENONIOO bound, when the scan is performed with the 
theoretical uncertainty r = 10 x cjpgo we label this with "X100(r=10)" ; conversely, when we apply 
r = 1 X (Tp9o we use the label "X100(r=l)". Similarly, for the dSphs bound, when the scan is 
performed with the theoretical uncertainty r = 50 x (o"f )tot,95 estimated for the choice of the DM 
profile we use "dSphs (Unc.)", whereas when we select the NEW profile we use "dSphs (NEW)". 

5.1 Impact on the CMSSM parameter space 

To start with, in Eig. la we show the marginalized posterior pdf in the {ruQ, mi/2) plane for the set 
of: non-LHC set of constraints, the LHC ax constraint from CMS analysis of 1.1 fb^^ data and the 
XENONIOO 90% C.L. bound on with a conservative estimate of the theory error r = 10 x fXp gQ. 
Our intention here is to summarize the present impact of experimental limits on the PS of the 
CMSSM before we move to including the ID observations and sensitivities. We show the 68% C.L. 
credible regions in dark blue, and the 95% C.L. regions in light blue. The black line indicates the 
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Figure 1: Marginalized posterior pdf in the {'iTLQ,mi/2) plane. The 68% C.L. credible regions are 
shown in dark blue and the 95% C.L. regions in light blue, (a) The CMSSM parameters constrained 
by the non-LHC, ax and XENONIOO data, (b) The same constraints as in (a), but the theoretical 
uncertainty on XENONIOO is strongly reduced (r = 1). The solid black line shows the 95% C.L. 
exclusion bound from the ax analysis of CMS. 



CMS 95% C.L. exclusion limit ^Ij, although we stress that, as previously in [9|, we apply it through 
our own approximate likelihood function, as discussed in Sec. [3) 

The relative impact of the constraints applied in Fig. la has been extensively analyzed in [9], 
for a smaller range of mo and 'mi/2- In particular one can see the two distinct 68% C.L. credible 
regions already found in [9j. One of these regions peaks in the stau-coannihilation/A-resonance 
region, where mo is relatively small and mi/2 ^ 500 GeV. (The ^-resonance actually extends to 
the much larger 2a region.) The other appears in the FP/HB region at large uiq and extends to 
much larger values of tuq than those considered in jSJ. As noted in [9j, the effect of the ay bound 
is to mostly push the high-probability credible regions up and outside the experimental exclusion 
countour. 



Figure la includes the 90% C.L. bound on cjp^ with a conservative theoretical uncertainty. 



10 X CTp 9Q, as explained in Sec. 



3.2 



We have checked that removing the bound from the 
figure would make hardly any difference. The near independence of Fig. la on the XENONIOO 
limit reflects the fact that at present the LHC limits are much stronger, while DD experiments are 
still marred by large theoretical uncertainties. We thus confirm and extend to larger parts of the 
FP/HB region the finding of Ref. [9] that the further effect of the current XENONIOO bound on 
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the posterior is minimal. 



In Fig. lb we investigate the effect of assuming a highly optimistic uncertainty of r = cr^ gg, as 



described in Sec. 3.2 We can see in Fig. lb that it strongly affects the FP/HB region by reducing 
the area of the 95% C.L. credible region above mi/2 ^ 200 GeV at large mo. The reason for this 
is clear: in the FP/HB region the higgsino component of the neutralino increases, and the SI cross 
section with nuclei is enhanced, it thus becomes more constrained by DD limits. On the other hand, 
it appears that even by reducing the theoretical uncertainties to only 100% of the experimental 
upper bound (a highly unlikely possibility) the impact on the PS is still not sizeable. Below we will 
perform a similar analysis in the (m^,(Tp^) plane. 
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Figure 2: Marginalized posterior pdf in the {mQ,mi/2) plane. The 68% C.L. credible regions are 
shown in dark blue and the 95% C.L. regions in light blue, (a) The CMSSM parameters constrained 
by the non-LHC, ax, XENONIOO and dSphs data. The uncertainty on the DM halo profile is taken 
into account, (b) The same constraints as in (a), but the DM halo profile is chosen to be the NFW 
one. The solid black line shows the aj- 95% C.L. exclusion bound. 

Next, in Fig. 2 we show the effects of the dSphs bound from the FermiLAT Collaboration on 
7-rays on the marginalized posterior pdf in the {mQ,mi^2) plane. As we explained in Sec. 



4.1 



we 



have decided to adopt in our analysis two different cases, similarly to what we did for DD. First, 
to be on a conservative side, we take into account the uncertainties due to the J-factor, related to 
a freedom in selecting the halo profiles. Secondly, we assume the J-factor to be due to the NFW 
DM profile, following the approach adopted in the experimental analysis |14j . 

Figure 2a shows the result of taking into account uncertainties related to the choice of the DM 
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profile. As expected, imposing a theoretical error of 50 times the given bound does not have much 
visible effect on PS, as can be clearly seen by comparing Fig. 2a with Fig. la, where the constraint 
from ID is absent. On the other hand in Fig. 2b we show the result of assuming the NFW profile 
for DM distribution. In this case, the bound from dSphs cuts significantly into the FP/HB region, 
and the impact is stronger than observed with XENONIOO alone, even for the idealized case of 
extremely low theoretical uncertainties on ap, depicted in Fig. lb. We can see that the bound from 
dSphs 7-ray fluxes, and the NFW profile assumed, significantly lowers the statistical relevance of 
the FP region at mo ~ 2500 GeV, mi/2 ~ 400 GeV. This is because, in this region, the enhanced 
higgsino component of the neutralino produces larger {crv), which for a cuspy NFW profile results 
in higher fluxes which in turn are more strongly constrained by the FermiLAT data. 
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Figure 3: The same as in Fig. [2] but in the (^0;tan/3) plane. 



The effect of the above constraints on the marginalized posterior pdf in the (Aq, tan/3) plane is 
presented in Fig. 3. First, we have checked that, without including the dSphs 7-ray constraint, our 
results agreed well with those of Ref . [9] . We direct the reader to that paper for a detailed discussion 
of the impact of these limits on the PS. Including the dSphs 7-ray constraint, but conservatively 
without selecting a specific DM profile (Fig. 3a), has a very limited effect on the posterior, similarly 
to what we have found for the {TnQ,mi/2) plane above. Hoewever, some differences with Fig. 8b 
of [9] can be seen. Since our scan spans a wider fraction of the FP/HB region, the region around 
tan/3 ~ 30-40 gains higher statistical significance than reported in that paper, and it is contained 
in a 95% C.L. credible region. This behavior is due to the much enlarged Icr-credible region at 
large niQ and mi/2 ~ 300 GeV that can be observed in Fig. la. This region does not show up 
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in [9] at the la level but only at the 2a level. As explained above, its presence here means that the 
experimental constraints now favor a larger fraction of the FP/HB region which features a higher 
higgsino content in the neutralino, and therefore a larger significance of tan /3 across a wide range 
of values 25 < tan/3 < 55 [9j, which enhances also the region around 30-40 mentioned above. 

On the other hand, we can see in Fig. 3b that, assuming the NFW profile for DM distribution 
in interpretting the data from dSphs 7-ray flux, has the effect of significantly constraining the 
scattering and annihilation cross sections, which reduces large fractions of high probability regions 
of the the FP region, and also pushes back the statistical significance of part of the region at 
tan/3 ~ 40 to less than 95% C.L. 

5.2 Effects on the observables 

We can better understand the effect of including the FP/HB region on our scan by looking at the 
posterior distribution of the observables relevant for DM detection: 0"^^ (erf) and the neutralino- 
proton SD cross section cTp^. 

In Fig. 4a we show the posterior pdf projected on the {my,,ap^) plane, with the non-LHC, ax 
and XENONIOO (assuming r = 10 x a^^gQ) constraints imposed. The solid black line shows the 
XENONIOO 90% C.L. bound, although we again emphasize that we implement it with a smearing 
through our likelihood function, and not as a hard cut. We can recognize three familiar high 
probability regions. A roughly horizontal one corresponds to the FP/HB region. One (mostly 
2o"-credible region but with a narrow la piece) decreasing with is given by the A-resonance 
effect. Finally, an "island" of high probability at 200GeV < < 400GeV corresponds to the 
stau-coannihilation region. 

A comparison with the corresponding figure in [9] shows the emergence of a sizeable Icr-credible 
region above the experimental bound, for lower < GeV. As stated above, this is a direct 
consequence of our extension to include a much larger part of the FP/HB region. Figure 4a 
confirms that the XENONIOO bound itself does not affect much the statistical significance of this 
region of PS, due to the large theoretical uncertainties associated with The idealized case 
where those uncertainties are reduced is shown in Fig. 4b where the 68% C.L. credible region above 
C7p gQ shrinks somewhat but is still present there. 

The effect of applying the dSphs bound is shown in Fig. 4c where we assume a large uncertainty 
in the DM profile, and in Fig. 4d assuming the NFW profile. Unsurprisingly, we can see that 
adopting conservative (and realistic) uncertainties on either ap^ or the dSphs 7-ray fiux has a 
rather limited effect beyond the LHC limit. On the other hand, the lower the error involved in 
defining the observable is assumed, the higher becomes the capability of the experimental limits 
to actually make an impact on the CMSSM PS. But more importantly, a comparison of Figs. 4b 
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Figure 4: Marginalized posterior pdf in the (m^jap^) plane. The 68% C.L. credible regions are 
shown in dark blue, the 95% C.L. regions in light blue, (a) The impact of the non-LHC, a-r 
and XENONIOO constraints, (b) The same constraints as in (a), but the theoretical uncertainty 
on XENONIOO is strongly reduced, (c) The impact of non-LHC, ar, XENONIOO and dSphs 
constraints. The uncertainty on the DM halo model is taken into account, (d) The impact of the 
non-LHC, ot, XENONIOO and dSphs constraints. The uncertainty on the DM halo model is not 
considered. The solid black line shows the XENONIOO 90% C.L. exclusion bound. 



and 4d shows that the effects of direct and indirect detection searches at present sensitivities would 
become comparable if the uncertainties involved had a reduced size. This fact was already hinted at 
when comparing Figs, lb and 2b, and shows how completely different experiments can in principle 
be competitive in testing the CMSSM. 

Figure 5 shows the marginalized posterior pdf in the {rriy., {cv)) plane. To start with, the effect 
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Figure 5: Marginalized posterior pdf in the (m^, (iTw)aim) plane. The 68% C.L. credible regions are shown 
in dark blue, the 95% C.L. regions in light blue, (a) The impact of non-LHC and ax constraints only, (b) 
The impact of non-LHC, ar and XENONIOO constraints. The theoretical uncertainty on XENONIOO is 
reduced, (c) The impact of non-LHC, ar and dSphs constraints. The DM halo follows the NFW profile, 
(d) The impact of non-LHC, ar, XENONIOO and dSphs constraints, with the NFW profile for DM halo 
assumed. 



of applying the non-LHC and ax constraints only is presented in Fig. 5a where we can see a familiar 
high probability horizontal region at ((Tv)ann ~ 3 x 10^^^ cm^'^ which arises in order to satisfy 
the relic abundance constraint in the FP/HB (and bulk) region from neutralino pair-annihilation. 
Additionally, there is an "island" of high probability at 200GeV < < 400GeV and much lower 
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(o'f)ann Corresponding to the stau-coannihilation region, where the neutraUno annihilation cross- 
section does not provide the only contribution to the relic density. 

Figure 5b shows the most optimistic scenarios where the XENONIOO bound is applied with 
the idealized reduced theoretical uncertainties and likewise Fig. 5c shows the impact of the dSphs 
bound with the cuspy NFW profile only, and without the XENONIOO applied. Finally, Fig. 5d 
shows the impact of dSphs when added to the non-LHC, a-r and XENONIOO constraints. It is clear 
that, when considered with the NFW profile, the dSphs bound disfavors the lu-credible region at 
large (av) around and just above mw, as this is the region that produces the highest 7-ray flux 
through the ratio {(Tv)a.nn/fn'^ in Eq. ([t]). On the other hand, the Icr-credible region at > mw 
is still unconstrained by DD and ID experiments, even with idealized low uncertainties assummed. 



5.3 Effects due to IceCube 
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Figure 6: Marginalized posterior pdf in the (m^,(Tp^) plane. The 68% C.L. credible regions are 
shown in dark blue, the 95% C.L. regions in light blue, (a) The impact of the non-LHC, a^, 
XENONIOO and dSphs constraints with the NFW profile, (b) The impact of non-LHC, a^, 
XENONIOO and dSphs constraints. The theoretical uncertainty on XENONIOO is reduced, and 
the DM halo follows the NFW profile. The solid line is the one-year 95% C.L. sensitivity in con- 
tained events at IceCube, the dashed line the corresponding sensitivity for DeepCore, and the 
dash-dotted line the sensitivity in upward events. 



In this subsection we analyze the effect on the CMSSM of the projected one-year 95% C.L. ex- 
perimental sensitivities for the 80-string configuration at IceCube, and the impact of the additional 
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six DeepCore strings. In the assumption of equilibrium in the Sun, the event rates depend directly 
on the neutralino capture rate, and consequently are most strongly dependent on the SD scattering 
cross section (although for some regions of PS the SI cross section can be important [53j)- 

In both panels of Fig. 6 we plot the marginalized pdf in the (m^,(Tp^) plane, after applying 
all the constraints discussed above. We only consider the NFW profile, which gives the more 
stringent bounds from FermiLAT. In Fig. 6a we show the case including the large uncertainties 
for DD searches, and in Fig. 6b the one with reduced uncertainties. The three black lines show 
the projected sensitivities of the IceCube detection channels considered in this work. The dash-dot 
line shows the sensitivity for upward events at IceCube, the solid line the sensitivity in contained 
events, and the dashed line the sensitivity for DeepCore. We confirm the complementarity of the 
three configurations to test different ranges in the neutralino mass |44| I45j . Due to its denser array, 
DeepCore has a lower energy threshold, so that it is better suited to probe the region at low mass, 
niy^ ~ 100 GeV. As the DM mass increases, though, the neutrino flux becomes suppressed by l/m^, 
so that the bigger volume available in the IceCube 80-string configuration becomes more important. 
The result is increased sensitivities in the larger detector, which can overcome the observational 
capabilities of DeepCore. 

In the first three panels of Fig. 7 we show the posterior probability in the {m^, A'events) plane, 
where the muon rates at the detector, N^^^^ts (Fig- 7a), A^'^^^^-g (Fig. 7b) and A^^^^-g (Fig. 7c) are 



given in Eqs. (12), (13) and (14), respectively. In dark blue we can spot three Ic-credible regions 
in the mass range 200GeV < niy. < 400GeV. Since they correspond to the two high-probabilty 
regions depicted in Fig. 6 for the same mass range, the splitting of the lower regions indicates some 
dependence on SI interactions. 

For contained events in IceCube, and for DeepCore, about 20-25% of the 2cr-favored PS is 
above the projected sensitivity, thus indicating that the next round of data from the IceCube 
Collaboration has the potential of significantly constrain the PS. In Fig. 7d we show the region of 
PS in the (mo, ^1-1/2) plane which is above the sensitivity at DeepCore. It overlaps with the region 
tested by the FermiLAT 7-ray flux from dSphs. 

Some very enticing possibilities arise in case IceCube/DeepCore lived up to its expected dis- 
covery potential. In particular, we found that if the projected sensitivities were echieved, and the 
experimental and theoretical uncertanties were kept reasonably small, in one year time DeepCore 
could rule out the fraction of the FP region which is mostly responsible for neutralino annihilation 
to W^W~ at the 95% C.L. To illustrate this point we show in Fig. 8 the impact of some of the exclu- 
sion bounds and projected sensitivies considered in this work on models with a single predominant 
neutralino annahilation channel. Each point represents a case for which neutralinos annihilate 90% 
of the time to the selected final state, and for which the bound on the relic abundance is respected 
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Figure 7: (a)-(c) Marginalized posterior pdf in the (m^^, A'events) plane. The 68% C.L. credible 
regions are shown in dark blue, the 95% C.L. regions in light blue, (d) Marginalized posterior 
pdf in the (771-0,^^/2) plane with DeepCore one-year 95% C.L. sensitivity (in red dashed lines). 
In all plots the constraints are non-LHC, ot, XENONIOO and dSphs. The DM halo follows the 
NFW profile. Upward events at IceCube are shown in (a), contained events in (b), and events at 
DeepCore in (c). 



(0.05 < < 0.20). The likelihood used to draw the point distribution in space is subject to all 

the constraints discussed above. The DM profile is assumed to be the NFW. We show in Fig. 8a 
the impact of the contained event sensitivity (black horizontal line) on these models. It is more 
constraining than the upward event sensitivity (not shown), but no final state can be exluded in one 
year time of observation. On the other hand. Fig. 8b shows the same analysis for DeepCore. We 
can see that the T^VF-channel could be excluded at the 95% C.L., and so could any CMSSM-like 
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Figure 8: The impact of some of the exchision bounds and projected sensitivities considered in 
this work on models with a predominant neutrahno annihilation channel, (a) Projected one-year 
95% C.L. sensitivity in contained events, (b) Projected sensitivity at DeepCore. (c) LHC ay 
bound, (d) XENONIOO bound. 

model which for which the DM annahilates into W~^W~ . The advantage of DeepCore comes from 
the lower energy threshold. Since the W^H^-channel opens up around myy, which is very close to 
the threshold predicted for the IceCube larger volume, there is not enough phase space over which 
the muon flux can be integrated in the contained events. Consequently the number of hits at the 
detector will not be sizeable. On the other hand, since the threshold at DeepCore is ~ 35 GeV, even 
neutralinos that annahilate into the M^VF-channel can give enough muon events to be succesfully 
detected. We can see in Fig. 8c that the FP region is mostly responsible for cases with predominant 
DM annihilation to W~^W~. We show in Fig. 8d the corresponding analysis for the XENONIOO 
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exclusion bound on (black line). A fraction of the cases that annahilate to W~^W are below 
the bound so nothing definitive can be said about the exclusion potential in any channel. 

6 Summary 

In this paper we investigated the impact of the present bounds and projected sensitivities from DM 
indirect detection searches on global Bayesian inferences of the CMSSM. We applied the 95% C.L. 
upper bound on (av) from the recent analysis of 7-ray fluxes from dSphs by the FermiLAT Col- 
laboration, and calculated the effect of the projected one-year 95% C.L. sensitivities for neutrino 
fluxes from WIMP annihilation in the Sun at IceCube and DeepCore. We also applied the con- 
straints used in Ref . [9] , with respect to which we extended our scan deep into the FP region, up to 
niQ = 4 TeV and m.1/2 = 2 TeV. In particular, we considered the recent lower bounds on sparticle 
masses from the CMS ut search with 1.1 fb^^ integrated luminosity, and the recent 90% C.L. upper 
bound on from XENONIOO. With respect to [9], we improved the code for determination of 
the efficiency maps, by simulating the CMS detector response with PGS4. 

We found that the impact of the bounds from dSphs is very sensitive to the choice of the DM 
halo profile, and that when the theoretical uncertainties are fully taken into consideration the con- 
straining power of the present data is limited. This applies equally to the current bound from 
XENONIOO. On the other hand, we found that a less conservative, in fact currently idealized, 
treatment of the theoretical uncertainties involved, obtained often by including the experimen- 
tal uncertainties only (as given by the experimental collaborations) shows significant constraining 
power on the FP region. Even more interestingly, all DM constraints, both from DD and ID, become 
competitive in this case. Therefore, we incite the experimental and theoretical groups to renew 
their efforts to abate those uncertainties. Finally we showed that the projected one year sensitivity 
at DeepCore has the potential to further constrain the FP region, and to exclude the fraction of 
PS over which DM annihilates predominantly to W~^W~ final states. Thus, the next round of data 
from IceCube/DeepCore has the potential to rule out any new physics models whose coupling to 
SM particles is comparable in strength to the CMSSM, and whose products of annihilation are 
M^-bosons. 
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